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1 Introduction 

The first proton run at a centre-of-mass energy (y/s) of 7 TeV at the Large Hadron Collider 
(LHC) took place from 30th March 2010 until 4th November 2010, with almost 50 ptr 1 of 
integrated luminosity delivered to the ATLAS and CMS detectors, half of it in one week 
at the end of October. The first measurement of W and Z cross sections was reported by 
both ATLAS and CMS at the ICHEP conference in July 2010, and subsequently published 
with integrated luminosities of approximately 320 nb -1 for ATLAS [1] and 2.9 pb _1 for 
CMS [2]. Preliminary results with the full 2010 data set were presented in March 2011 
by both ATLAS [3] and CMS [4]. (The LHCb experiment has also presented preliminary 
results for the W and Z cross sections in the forward region [5].) The observation of 
the first few top candidate events was also reported at the ICHEP conference, and first 
ti cross-section measurements were published with integrated luminosities of 2.9 pb _1 for 
ATLAS [6] and 3.1 pb" 1 for CMS [7]. Again, preliminary results with the full 2010 data 
set were presented in March 2011 [8, 9] and CMS have since published their measurements 
in the dilepton channel [10] and the lepton+jets channel [11]. The precision of the cross 
sections for W, Z and ti production will improve substantially in the near future as the 
2011 LHC run progresses at an impressive pace, with the remarkable milestone of 1 fb~ x 
of integrated luminosity already delivered to the ATLAS and CMS detectors by 17th June 
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2011. Searching for, or excluding, the Standard Model Higgs boson (H), produced mainly 
from gluon-gluon fusion through a top-quark loop, requires precise knowledge of the the- 
oretical cross section. All these benchmark processes (W, Z, ti and gg — > H) are sensitive 
to the parton distribution functions (PDFs) of the proton. The proton PDFs are deter- 
mined by several groups from (global) analysis of a wide range of deep-inelastic scattering 
(DIS) and related hard-scattering data, taken at the HERA and Tevatron colliders and 
from fixed-target experiments at lower centre-of-mass energies; see refs. [12-14] for recent 
reviews. 

The aim of this paper is to compare next-to-leading order (NLO) and next-to-next-to- 
leading order (NNLO) predictions, calculated within the standard framework of collinear 
factorisation, for the benchmark processes using the most recent publicly available PDFs 
from each group, all available using the current lhapdf V5.8.5 [15] interface. We will pay 
particular attention to the values of the strong coupling, as, used by each fitting group 
and to the corresponding as dependence of the benchmark cross sections. We do not 
aim to come up with a single "best" prediction for each of the benchmark cross sections 
together with a complete evaluation of all sources of theoretical uncertainty. We consider 
only uncertainties due to PDFs and as on fixed-order (NLO and NNLO) predictions. We 
do not consider, for example, optimal (factorisation and renormalisation) scale choices and 
variations, electroweak corrections, or the effect of threshold resummation. For the gg — )■ H 
process, we do not consider (Ca tt as) "-enhanced terms, use of a finite top-quark mass in the 
calculation of higher-order corrections, bottom-quark loop contributions, etc. Discussion 
of these issues can be found in the recent Handbook of LHC Higgs Cross Sections [16]. 
The findings regarding PDFs and as reported here will be relevant also for more complete 
calculations, and the PDF and as uncertainties often make up a sizeable or dominant part 
of the total theoretical uncertainty. We will present results only for total cross sections 
rather than differential distributions, and we do not attempt to account for the precise 
experimental cuts, although for fairly inclusive cuts the results should be similar. The 
findings for y/s = 7 TeV will be similar also for slightly higher energies, should this happen 
before the planned LHC shutdown in 2013. Significant updates to the PDFs are expected 
by the time the LHC runs at up to \/s = 14 TeV, expected from 2014, therefore we do not 
present results at that centre-of-mass energy in this paper. 

The background leading up to this paper is as follows. The LHC Higgs Cross Section 
Working Group [17] was formed in January 2010 and requested a recommendation from 
the PDF4LHC Working Group [18] on PDFs and as values (and their uncertainties) to 
be used for cross-section calculations. A proposal was made, initially by J. Stirling, to use 
different PDFs to calculate LHC benchmark processes at 7 TeV, specifically production 
of W ± , Z°, ti and gg -)• H for M H = 120,180,240 GeV. The aims were to establish 
the degree of compatibility and identify outliers amongst PDF sets, and to compare cross 
sections at the same as values, thereby showing to what extent differences in predictions 
are due to the different as values adopted by each group, rather than differences in the 
PDFs themselves. The results at NLO, initially presented by G.W. in a PDF4LHC meeting 
at CERN on 26th March 2010 [19], formed the basis for the subsequent PDF4LHC Interim 
Report [13] and PDF4LHC Interim Recommendations [20] used in the Handbook of LHC 
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Table 1. Comparison of major PDF sets considered, and their gross features distinguished by 
the main classes of data included (upper part of table) and important aspects of the theoretical 
treatment (lower part of table), specifically regarding the treatment of heavy quarks in DIS and the 
provision of NNLO PDFs. More refined differences between PDF sets are described in the text. 

Higgs Cross Sections [16]. 

In this paper, we extend 1 the comparison of benchmark cross sections in a number of 
ways compared to the original PDF4LHC exercise [13, 19]. We will account for new NLO 
PDF sets released since March 2010, particularly the new CT10 [24] and NNPDF2.1 [25] 
analyses, and we extend all comparisons to NNLO, 2 making use of the recently released 
hathor public code [27] for an approximate NNLO ti calculation. A companion web- 
page [28] will be used as a repository for supplementary plots and may be updated after 
the publication of this paper, for example, as new PDF sets continue to be released. It is 
hoped that this paper and the companion webpage will prove to be a useful resource to 
accompany the publication of precise measurements of LHC cross sections from the 7 TeV 
proton runs. In a separate but related paper [29], we comment on the ability of the Teva- 
tron jet data to discriminate between PDF sets with different high-x gluon distributions 
and we respond to some recent criticism [26, 30-34] regarding theoretical uncertainties due 
to PDFs (and as) in the gg — >• H cross section at the Tevatron (and LHC). 

The structure of this paper is as follows. In section 2 we briefly describe the PDF sets 
considered in the study and specify formulae used for uncertainty calculations, in section 3 
we compare partonic luminosities, in section 4 we plot total cross sections versus as(M^), 
in section 5 we make some comments on the PDF4LHC recommendations, in section 6 we 
discuss W and Z production 3 in more detail, then finally we conclude in section 7. 

2 PDF sets considered 

We will consider only public sets, defined to be those available for use with the latest 
lhapdf V5.8.5 [15]. We consider NLO PDFs from MSTW08 [37], CTEQ6.6 [38], CT10 [24], 
NNPDF2.1 [25], HERAPDF1.0 [39], HERAPDF1.5 (preliminary) [40], ABKM09 [41] and 
GJR08 [42, 43]. We consider NNLO PDFs from MSTW08 [37], HERAPDF1.0 [39], ABKM09 [41] 
and JR09 [44, 45]. The broad distinctions between data sets fitted and aspects of the theo- 
retical treatment are summarised in table 1; see also the descriptions in refs. [12-14]. Only 
three groups (MSTW08, CTEQ6.6/CT10 and NNPDF2.1) make fully global fits to HERA 
and fixed-target DIS data, fixed-target Drell-Yan production, and Tevatron data on W, Z 

1 Much of the new material was initially presented in some recent talks [21-23]. 

2 See also ref. [26] for NNLO comparisons of W, Z and Higgs production with different PDFs. 

3 The PDF dependence of W and Z production at the LHC was also previously discussed in refs. [35, 36]. 
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and jet production, although GJR08 includes all these processes other than Tevatron W 
and Z production. The HERAPDF1.0 fit includes only the combined HERA I inclusive 
data, while the HERAPDF1.5 fit additionally includes the preliminary combined HERA 
II inclusive data. The CT10 and NNPDF2.1 global fits include the combined HERA I 
inclusive data [39], while the other global fits include the older separate data from HI 
and ZEUS. The MSTW08, CT10, NNPDF2.1 and GJR08 fits include Tevatron Run II 
data, while CTEQ6.6 uses only Tevatron Run I data. Most groups now treat the heavy- 
quark contribution to DIS structure functions using a general-mass variable flavour number 
scheme (GM-VFNS), other than ABKM09 and GJR08/JR09 who use a fixed flavour num- 
ber scheme (FFNS). The change from the inadequate zero-mass variable flavour number 
scheme (ZM-VFNS) to the GM-VFNS was the major improvement between NNPDF2.0 [46] 
and NNPDF2.1 [25], now allowing a meaningful comparison to other NLO global fits. The 
NNPDF fits parameterise the starting distributions at Qq = 2 GeV 2 as neural networks and 
use Monte Carlo methods for experimental error propagation. The other groups all use the 
more traditional approach of parameterising the input PDFs as some functional form in x, 
each with a handful of free parameters, and use the Hessian method for experimental error 
propagation with differing values of the tolerance A% 2 , that is, the change in the global 
goodness-of-fit measure relative to the best-fit value. Contrary to the "standard" input 
parameterisation at Qq > 1 GeV 2 , the GJR08/JR09 sets use a "dynamical" input param- 
eterisation of valence-like input distributions at an optimally chosen Qq < 1 GeV 2 , which 
gives a slightly worse fit quality and lower as values than the corresponding "standard" pa- 
rameterisation, but is nevertheless favoured by the GJR08/JR09 authors. Public NNLO fits 
are only available from MSTW08, HERAPDF1.0, ABKM09 and JR09. (The first NNLO 
fits from the CTEQ/CT and NNPDF groups should be available soon, together with HER- 
APDF1.5 NNLO, but these analyses are not yet published.) The Tevatron jet cross sections 
are excluded from the JR09 fit, where complete NNLO corrections are unavailable, whereas 
they are included in the MSTW08 NNLO fit by making the approximation of using the 
NLO partonic cross section supplemented by 2-loop threshold corrections [47]. 

2.1 Definition of PDF uncertainties 

All groups other than NNPDF propagate experimental errors on the fitted data points by 
diagonalising the nxn covariance (or Hessian) matrix, with n free parameters, and provide 
2n eigenvector PDF sets, (k = 1, . . . , n); see the discussion in ref. [13] for more details. 
We use asymmetric uncertainties (apart from for ABKM09, where St = S^) to calculate 
the PDF uncertainty, denoted (AF)±, on an observable F such as a hadronic cross section: 



(AF)_ 



{max [ F(S+) - F(S ), F(S~) - F(S ), 0] } 2 , (2.1) 



\ k=l 



(AF)_ = Yl i max [ F ( S °) ~ F ( S k)> F ( S °) ~ F ( S k)> 0] } 2 > (2-2) 
\ k=i 

where n = 20, 22, 26, 10, 25, 13 for MSTW08, CTEQ6.6, CT10, HERAPDF1.0/1.5, ABKM09, 
GJR08/JR09, respectively. The HERAPDF1.0/1.5 fit provides additional "model" and 
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"parameterisation" asymmetric errors, which are added in quadrature to the "experimen- 
tal" errors. The HERAPDF1.0 NNLO fit does not include any uncertainties, only the 
central fit for two values of as(M§) = 0.1145 and 0.1176. 

For NNPDF2.1, instead the average (F) and standard deviation AF are taken over 
the observable F calculated with each PDF replica set, Sk (k = 1, . . . , N rcp = 100), that is, 

N lep 

rep k=i 



AF = \/^T (<F2) - <F>2) - (2 ' 4) 

For MSTW08, eigenvector PDF sets are consistently provided for both the estimated 
68% confidence-level (C.L.) and 90% C.L. uncertainties. For other groups, where usu- 
ally only 68% is provided (but only 90% for CTEQ6.6/CT10) we assume that the 68% 
C.L. and 90% C.L. uncertainties are related by a factor of 1.64485. For NNPDF2.1, while 
it is possible to calculate exact confidence levels (within the limited statistics), we will 
assume that a standard deviation, AF given by eq. (2.4), represents a 68% C.L. uncer- 
tainty and then rescale by 1.64485 to get a 90% C.L. uncertainty where necessary. We 
will sometimes only show 68% C.L. plots in this paper, but the corresponding plots at 90% 
C.L. can be seen at ref. [28]. The PDF uncertainties from HERAPDF1.0/1.5 and ABKM09 
include contributions from variation of the heavy-quark masses, m c and mb, while these 
are not included for CTEQ6.6, CT10 and GJR08/JR09. For MSTW08 and NNPDF2.1, 
the uncertainty from the choice of m c and will not be included here, although it can be 
included in principle according to the prescriptions presented in refs. [25, 48]. Variations of 
m c = 1.40 ± 0.15 GeV led to a change of just over 1% in the W and Z total cross sections 
for MSTW08, while varying m& = 4.75 ± 0.25 GeV gave a negligible change (0.1%), such 
that the m c ^ uncertainties are not important when added in quadrature with the 2-3% 
"PDF+as" uncertainty [48]; see also tables 3 and 4 in section 6. Similar variations of m c 
and mi, for Standard Model Higgs boson production (Mg = 120 GeV) at the 7 TeV LHC 
had essentially no impact [48]. Additional uncertainties from the particular GM-VFNS 
variant are estimated [49] to be ~ 2% at NLO, but less than 1% at NNLO, and will not 
be considered here. Combined HERA I data [39], not included in the MSTW08 fit, gave 
an upwards shift in the W and Z total cross sections by 1-2%, with a smaller change 
in the Higgs cross section [50]. The combined HERA I and II inclusive data, together 
with the combined charm structure function _pc harm 5 j-^ included in a future update to 
MSTW08, after it is finally published. 

2.2 Definition of PDF correlations 

Defining a correlation cosine between two quantities F and G in the Hessian approach, 

1 " 

COS ^ FG = 4AFAG £ " F{S ^ [G( ^ +) " ° {S ^ ' (2 - 5) 

k=l 
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where the uncertainties AF and AG are calculated using the symmetric formula: 



AF - - 

2 



£[F(S+)-F(S t 7)] 2 , (2.6) 

k=l 



then values of cos <j>FG ~ 1 mean that F and G are correlated, values of ~ — 1 mean that 
they are anticorrelated, while values of ~ mean that they are uncorrelated. In the Monte 
Carlo approach used by the NNPDF group, the equivalent of eq. (2.5) is given by 

, jV rcp (FG)-(F)(G) 

C0S4>FG= N TCP -1 AFAG ' (2 ' 7) 

where the average (F) and standard deviation AF are defined in eqs. (2.3) and (2.4), 
respectively. In both the Hessian and Monte Carlo approaches, an ellipse in the F—G plane 
can then be defined by the two parametric equations: 

F = F{Sq) + AF cos 9, (2.8) 
G = G{S ) + AG cos (6 + 4>fg) , (2.9) 

where 9 € [0, 27r] and F(Sq) should be replaced by (F) in the Monte Carlo approach. We 
give explicit examples in section 6 for the W and Z total cross sections at the LHC. Finally, 
we note that the uncertainty in a general function of the two quantities, H(F, G), is given 
by (see also section 4 of ref. [13]): 

,dH\ 2 ( A „dH\ 2 . , dHdH 



<A*> ={ AF W) + { AG 8g) +2AFAGW fG __. (2.10) 

As a specific example, the fractional uncertainty in the ratio of the two quantities R = F/G 
is therefore given by 

AZA 2 / AF\ 2 ( AG\ 2 AFAG 

so the uncertainty in the ratio can be obtained from the correlation cosine (or vice versa). 
Of course, eqs. (2.8), (2.9), (2.10) and (2.11) hold much more generally than just for PDF 
uncertainties and can be applied to define an ellipse also for experimental data; see ref. [23] . 

2.3 Values and uncertainties of strong coupling as 

The values of as{Mg), and the corresponding uncertainties, for MSTW08, ABKM09 and 
GJR08/JR09 arc obtained from a simultaneous fit with the PDF parameters. Other groups 
choose a fixed value, generally close to the world average [51], and for those groups we 
assume an uncertainty of ±0.0012 (±0.0020) for a 68% (90%) C.L. [13]. The central values 
and l-o" uncertainties are depicted in figure 1 as the larger symbols and error bars, while 
the smaller symbols indicate the PDF sets with alternative values of as{M^) provided by 
each group. The fitted NLO as(M^) value is always larger than the corresponding NNLO 
as(Mg) value in an attempt by the fit to mimic the missing higher-order corrections, 
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(a) 



(b) 



NLO c/ s (M 2 ) values used by different PDF groups 
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■ CTE06.6 

▲ CT10 
O NNPDF2.1 
□ HERAPDF1 .0 
ABKM09 
GJR08 



a s (M*) 



NNLO a s (M ) values used by different PDF groups 



66% C.L. PDF 
£ M ST WO 8 

■ HERAPDF1 .0 
A ABKMQ9 

GJR08JR09 



Closed symbols: NNLO 
Open symbols: NLO 



<x s < M z) 



Figure 1. Values of as(M^), and their 1-er uncertainties, used by different PDF fitting groups 
at (a) NLO and (b) NNLO. The smaller symbols indicate the PDF sets with alternative values of 
«s(Mf ) provided by each fitting group. The shaded band indicates the world average as(M§) [51]. 



which are generally positive. Note that most of the fitted NLO as(M^) values (other than 
GJR08) agree with the world average [51], while MSTW08 is the only set in agreement at 
NNLO. 

The prescription for "PDF+as" uncertainties in the MSTW08 case was given in 
ref. [52]. We compare this prescription with the alternative prescription of adding in 
quadrature in table 2 for all the benchmark processes. (The settings for the cross-section 
calculations are given later in section 4.) Adding in quadrature for MSTW08 was proposed 
as a simplification in refs. [16, 20] and expected to give slightly smaller uncertainties. 
For CTEQ6.6/CT10 and HERAPDF1.0/1.5 we simply add as uncertainties in quadrature 
with the PDF uncertainties, as recommended by the fit authors. For NNPDF2.1, we use 
the prescription given in refs. [53, 54], that is, we take the number of replicas for each 

as(Mg) value from a Gaussian distribution with mean of 0.1190 and standard deviation 

(i) 

0.0012. Rounding to the nearest integer values, this gives iVrcp = 4,25,71,100,71,25,4, 
for a^ } (M|) = 0.116,0.117,0.118,0.119,0.120,0.121,0.122, then we take the averages and 

standard deviations over jV rep = J2j N ^l = 300 replicas. For each a^ } (M|) value, the 
(i) 

A^rcp replicas are chosen randomly from the full 100 replicas available. To eliminate depen- 

(i) 

dence on the particular subset of N^p replicas chosen, the averages and standard deviations 
over the N rep = 300 replicas are themselves calculated many (10 5 ) times, then the results 
averaged. This final averaging is supplementary to the official NNPDF prescription and 
can be considered as an improvement at little computational cost. For ABKM09 and 
GJR08/JR09, as(M^) is one of the fit parameters used to produce the eigenvector sets, 
which implies mixing of as(M^) with PDF parameters, so a "PDF+ag" uncertainty is ob- 
tained automatically and it is not straightforward to calculate a "PDF only" uncertainty. 
Care must then be taken that a different as value is used in the partonic cross section for 
each eigenvector PDF set. 
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+4.3% 
-3.9% 
+4.4% 
-4.2% 


+3.5% 
-4.3% 
+4.2% 
-4.1% 


+5.3% 
-5.4% 
+5.2% 
-5.4% 


+8.8% 
-8.3% 
+9.5% 
-9.0% 


+7.9% 
-9.0% 
+9.8% 
-9.1% 



Table 2. NLO and NNLO predictions for total cross sections using MSTW 2008 PDFs [37, 52] 
comparing PDF-only uncertainties and PDF+ag uncertainties computed using either the exact 
prescription [52] or simply adding as uncertainties in quadrature with the PDF-only uncertainties. 



3 Comparison of partonic luminosities 

To understand properties of hadronic cross sections, such as PDF uncertainties or the de- 
pendence on collider energy, it is useful to consider the relevant parton-parton luminosities 
(see, for example, refs. [55-60]). We define the dimensionless variable, r = s/s, where 
r G (0, 1] and \/~§ is the partonic centre-of-mass energy. A PDF of flavour a = q,q,g is 
denoted as f a (x,/j, 2 ), with x the momentum fraction and \x the factorisation scale. We will 
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consider the following combinations of partonic luminosity: 



ds 



1 




1 dx 



(3.1) 



s 



q=d,u,s,c,b 



ds 

9£gg 
ds 



s 



s 



I 




where qq is relevant for W and Z production, is relevant for Higgs and ti production, 
and GG is relevant for dijet production [56, 61]. Note that the dependence of dC/dr only 
on the scaling variable r = s/s rather than separately on s and s is (mildly) broken by 
the choice of factorisation scale u 2 = s. This means, for example, that ratios of partonic 
luminosities at \fs = 7 TeV plotted as a function of \fr will not be identical to the same 
ratios at y/s = 14 TeV. The logarithmic integration over x in eqs. (3.1), (3.2) and (3.3) can 
be interpreted as an integration over the rapidity y = ln(x/\/r) of the final state. We will 
only consider rapidity-integrated partonic luminosities relevant to the total cross section in 
this paper, but it would be straightforward to also study the rapidity-dependent quantities. 
Note that the qq luminosity does not specifically include the correct flavour combinations 
for or Z° production. More detailed studies would, for example, include the correct 
couplings of the vector bosons to quarks and antiquarks, or consider the specific ud (for 
W + ) or du (for W~) partonic luminosities. Nevertheless, the generic qq luminosity defined 
in eq. (3.1) should exhibit the general features of the more precise partonic luminosities 
relevant for the W and Z total cross sections. 

In figure 2 we show the NLO parton-parton luminosities as the ratio with respect to the 
MSTW 2008 NLO luminosities, for the LHC at \fs = 7 TeV. The inner uncertainty bands 
(dashed lines) for HERAPDF1.0 correspond to the (asymmetric) experimental errors, com- 
puted using eqs. (2.1) and (2.2), while the outer uncertainty bands (shaded regions) also 
include the model and parameterisation errors, including uncertainties on heavy-quark 
masses but not on as- We use the default as values for each set, shown in figure 1. 
Recall that it is not possible to separate the "PDF only" uncertainty for ABKM09 and 
GJR08/JR09, therefore the uncertainty bands for those sets also include the as uncer- 
tainty, and the uncertainty bands for ABKM09 also include uncertainties on heavy-quark 
masses. This is undesirable but unavoidable given that these groups do not provide PDF 
sets for fixed as (and fixed m Ci & for ABKM09). These additional sources of uncertainty 
may be quite significant; see, for example, refs. [48, 52] for the as and m c ^ dependence of 
the MSTW08 PDFs. The relevant values of \fs~ = M\y,z are indicated on the plots of the 
qq luminosities in figure 2(a,b). There is fairly good agreement for the global fits, but more 
variation for the other sets, particularly when moving away from the central ~ My/^z- 
For the gg luminosities in figure 2(c,d), the relevant values of = MH,2mt are indi- 
cated, where = 2rrit (with wit = 171.3 GeV) is the threshold for ti production. Again, 
there is reasonable agreement for the global fits, but more variation for the other sets, 
particularly at large s, where HERAPDF1.0 and ABKM09 have much softer high-x gluon 
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Figure 2. NLO parton-parton luminosities as the ratio with respect to MSTW 2008. 



distributions. The somewhat larger gg luminosity of CTEQ6.6/CT10 at large s compared 
to MSTW08 may be due to the inclusion of Tevatron Run I jet data, which are known 
to prefer a larger high-x gluon distribution than the more reliable Run II data [37], while 
the larger gg luminosity at small s may be due to the positive-definite input gluon con- 
straint of CTEQ6.6/CT10, which is not imposed in the MSTW08 or NNPDF2.1 fits. By 
the momentum sum rule, the CTEQ6.6/CT10 gg luminosity at the intermediate s values 
relevant for gg — > H and ti production is then slightly smaller than that of MSTW08 and 
NNPDF2.1. The GG luminosities in figure 2(e,f), relevant for dijet production, interpolate 
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Figure 3. NLO parton-parton luminosities from HERAPDF1.0/1.5 as the ratio to MSTW 2008. 



between being ^-dominated at small s and (/(/-dominated at large s. Plots comparing 
the qq and gg luminosities from MSTW 2008 with HERAPDF1.0 and HERAPDF1.5 are 
shown in figure 3. The qq luminosity from HERAPDF1.5 is somewhat closer to MSTW08 
at large s than HERAPDF1.0, while there is little difference in the gg luminosity in going 
from HERAPDF1.0 to HERAPDF1.5. 

In figure 4(a,c,e) we show the NNLO parton-parton luminosities as the ratio with 
respect to the corresponding MSTW 2008 NNLO luminosities. The two HERAPDF1.0 
NNLO curves shown are for both asiM^) = 0.1145 and 0.1176, where the latter value 
gives the smaller gg and GG luminosity, and the larger qq luminosity, at low s values. The 
NNLO trend between groups is similar to at NLO. In figure 4(b,d,f) we also compare the 
MSTW08 NNLO parton-parton luminosities to the corresponding MSTW08 NLO lumi- 
nosities. The qq luminosity differs between NLO and NNLO by an amount comparable to 
the experimental uncertainty, while there is less difference in the gg and GG luminosities, 
partially due to the larger as(M^) value used at NLO compensating for the missing NNLO 
corrections to the evolution; see figure 1(b). 

In figure 5 we show the fractional uncertainties of the NLO parton-parton luminosities. 
For the global fits (MSTW08, CTEQ6.6/CT10, NNPDF2.1) there is good agreement be- 
tween the fractional uncertainties from different groups, particularly at intermediate values 
of s. This is quite remarkable considering the very different methods used for parameterisa- 
tion and experimental error propagation, especially comparing the novel NNPDF approach 
to the more traditional methods used by the other groups. In figure 6 we show the frac- 
tional uncertainties of the NNLO parton-parton luminosities, again comparing MSTW08 
NNLO to (a,c,e) other NNLO sets and to (b,d,f) MSTW08 NLO. The NNLO fractional 
uncertainties are similar to those at NLO. Uncertainties from ABKM09 and GJR08/JR09 
are generally larger than for the global fits, but this may partly be due to the additional 
sources of uncertainty from as (GJR08/JR09) and both as and m c ^ (ABKM09). 
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Figure 4. NNLO parton-parton luminosities as the ratio with respect to MSTW 2008. 



4 Total cross sections as a function of as(M^) 

The cross-section calculations are all performed in a 5-flavour ZM-VFNS using correspond- 
ing PDFs evolved with a maximum of 5 flavours, 4 as used also for the parton-parton 

4 The default PDF grids provided by NNPDF2.1 are evolved with 6 flavours for Q 2 > m 2 . Moreover, the 
nearest grid points to rrit = 175 GeV are at Q = 156 GeV (n/ = 5) and Q = 187 GeV (n/ =6), meaning 
that interpolation between these two grid points does not have a well-defined rif value. Although a variant 
of NNPDF2.1 with the default as {M%) = 0.119 is now provided with 5-flavour evolution, grids for other 
as values are only provided with 6-flavour evolution. Therefore we use only the 6-flavour NNPDF2.1 grids. 
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Figure 5. Fractional uncertainties of the NLO parton-parton luminosities. 



luminosity plots in section 3. We want to isolate the PDF (and as) dependence, there- 
fore we use precisely the same code for all PDF sets with common settings. The PDF+as 
uncertainties are computed using the recommended prescription of each fitting group, sum- 
marised in section 2.3. No attempt is made to evaluate other theoretical uncertainties and 
a single scale choice (//r = up = M\y, Mz, MH,m t ) is adopted for each process. We only 
present total cross sections, not differential distributions. For the PDF4LHC benchmark 
exercise at NLO [13], a modified mcfm V5.7 package [62, 63] with agreed parameters and 
input files was produced by J. Campbell. Instead we use private code (obtained from 



-13- 



(a) 



£„(qq) luminosity at LHC fjs = 7 TeV) 



-~ 0.2 
_l 

O 0.15 

S 0.1 
S'o.os 

'm 



-0.05 
-0.1 
-0.15 



r 



-0.2 



MSTW08 NNLO 

HERAPDF1.0 

ABKM09 

JR09 



lifflteiBffi M'D Tirirr; 

— ~~T~r — 



W Z 



(b) 



E q (qq) luminosity at LHC ^ = 7 TeV) 



'10- 3 



10- 1 




- MSTW08 NLO 

- MSTW08 NNLO 



w z 




(c) 



gg luminosity at LHC frjs = 7 TeV) 



-~ 0.2 
_i 

0.15 

sp 

§ 0.1 

I'o.os 

ra 

t 
0) 

u 

1 -0.05 
1 -0.1 

o 

'■5-0.15 



-0.2 
10 



MSTW08 NNLO 

HERAPDF1.0 

ABKM09 

JR09 




120 180 240 l . 

10" 2 M H (GeV) tt 10" 1 




(e) 



vs / s 



GG luminosity at LHC (sis = 7 TeV), G = g + 4/9 £ g (q+t[) 
0.2 



0.05 


-0.05 
-0.1 
-0.15 



-0.2 . 
10' 





1 1 1 1 1 , , , , 






MSTW08 NNLO 






HERAPDF1.0 




t 

: \ s 


ABKM09 
JR09 


^ilfllli 














^LIjjjjj-ujj- 1 


^^^^^ 




i i i i I . ... 





10- 1 



(d) 



(f) 



gg luminosity at LHC (Js = 7 TeV) 




10" 2 M„(GeV) tt 10" 1 



GG luminosity at LHC (Vs = 7 TeV), G a g + 4/9 L (q+q) 

»T 0.2 r 1 1 " i i J 




vs/ s 



vs/ s 



Figure 6. Fractional uncertainties of the NNLO parton-parton luminosities. 



J. Stirling) for total cross sections at NLO, also used for the majority of plots in ref. [13], 
which is much faster and was checked against the official mcfm package. Although fairly 
precise total cross sections for W^, Z Q and ti production have now been measured at 
the LHC, we resist the temptation to directly compare with data, since in many cases a 
more extensive evaluation of other theoretical uncertainties is necessary, and most of the 
cross-section measurements with the full 2010 data set are still preliminary [3, 4, 8, 9]. 
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Figure 7. W ± and Z° total cross sections, plotted as a function of as(M§), at NLO. 



4.1 and Z° production 

We calculate and cross sections, multiplied by the appropriate leptonic branching 
ratio, in the zero-width approximation with no virtual photon (7*) included in the Z° cross 
sections, as in section 15 of ref. [37]. We use the PDG 2008 [64] electroweak parameters, 
with the CKM matrix taken from a global fit in the Standard Model with unitarity con- 
straints imposed, i.e. the central values of the CKM matrix in eq. (11.27) of ref. [64]. The 
leptonic branching ratios are taken to be B(W -> Iv) = 0.1080 and B(Z -> it) = 0.33658. 
The W and Z masses are taken to be M w = 80.3980 GeV and M z = 91.1876 GeV, the effec- 
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Figure 8. W ± and Z° total cross sections, plotted as a function of as(M^), at NNLO. 



tive weak-mixing angle sm 2 (6\y) = 0.23149 (defining the vector coupling of the Z boson to 
quarks), the Fermi coupling constant Gf = 1.16637 x 10~ 5 GeV -2 and the conversion con- 
stant 1 GeV~ 2 = 389379.304 nb. The NLO total cross sections (with MSTW08 PDFs) cal- 
culated using the private code (originally from J. Stirling) were checked against the MCFM 
package provided for the PDF4LHC benchmark exercise [62, 63] and agreement was found 
to within the numerical integration error (of at most 0.1%) for all processes apart from Z°. 
The discrepancy in Z° production was traced to an ambiguity in the benchmark specifi- 
cation, where the provided mcfm version [63] used sin 2 (6w) = 1 - (M w /M z ) 2 = 0.22265, 
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Figure 9. W + and W total cross sections, plotted as a function of as(M§), at NLO. 



whereas we used the effective angle sin 2 (9\y) = 0.23149 [64] which absorbs electroweak ra- 
diative corrections, leading to a Z° total cross section 1.4% lower than the provided MCFM 
version [63]. The NNLO corrections were calculated in ref. [65] (see also ref. [66] for a 
correction). 

We show the (= W + + W~) and Z° total cross sections (and their ratio), plotted 
as a function of as(M|), in figures 7 and 8, with (a,c,e) 68% C.L. uncertainties on the left 
and (b,d,f) 90% C.L. uncertainties on the right. We show the separate W + and W~ total 
cross sections (and their ratio), plotted as a function of as(M^), in figures 9 and 10. 



-17- 



(a) 



NNLO W* -> Tv at the LHC fc/s = 7 TeV) 




Vertical ferror bars 
Inner: PDF only 
Outer: PDF+u 



Closed symbols: NNLO 
Open symbols: NLO 



(c) 



« S (M|) 



NNLO W" -> l"v at the LHC (s/s = 7 TeV) 



J2 



t 

m 




Closed symbols: NNLO 
Open symbols: NLO 



(e) 



a s (M 2 z ) 



NNLO W*/W ratio at the LHC (vis = 7 TeV) 









1.5 






3 


1.48 


b 


in 

+i 


1.46 


DC 






1.44 




1.42 




1.4 




1.38 




Vertical error b; 
Inner: PDF only 
Outer: PDF+u 





68% C.L. PDF 




£ MSTWOS 
■ HERAPDF1 .0 




^ ABKM09 

GJR08/JR09 



Closed symbols: NNLO 
Open symbols: NLO 



0.13 

« S (M|) 



(b) 



NNLO W* -> Tv at the LHC (s/s = 7 TeV) 




Vertical ferror bars 
Inner: PDF only 
Outer: PDF+u„ 



Closed symbols: NNLO 
Open symbols: NLO 



0.13 
|2\ 



(d) 



c 



■1 

t 

m 



« s ( M z) 



NNLO W" -> l"v at the LHC £ls = 7 TeV) 




(D 



NNLO W*/W ratio at the LHC (Js = 7 TeV) 




Figure 10. VK + and total cross sections, plotted as a function of as(M|), at NNLO. 



The common format of all the plots in figures 7, 8, 9 and 10 is that the markers are 
centred on the default as(Mg) value and the corresponding predicted cross section of each 
PDF fitting group. The horizontal error bars span the as(M|) uncertainty, the inner 
vertical error bars span the "PDF only" uncertainty where possible (i.e. not for ABKM09 
or GJR08/JR09), and the outer vertical error bars span the "PDF+as" uncertainty The 
effect of the additional as uncertainty is small for W and Z production. The dashed lines 
at NLO or the solid lines at NNLO interpolate the cross-section predictions calculated with 
the alternative PDF sets with different as(M^) values provided by each group, represented 
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by the smaller symbols in figure 1. The NNLO plots in figures 8 and 10 also show the 
NLO predictions (open symbols and dashed lines) together with the corresponding NNLO 
predictions (closed symbols and solid lines) to explicitly demonstrate how the size of the 
NNLO corrections depends on both the ag(M|) choice and the PDF choice. The NNLO 
corrections are reduced by taking different a,s(M§) values at different perturbative orders. 

The global fits are in good agreement for and Z° production, as was apparent from 
the qq luminosity plots in section 3, and there are small uncertainties in predictions for 
the W^/Z ratio. The W^/Z° ratio is insensitive to NNLO corrections (and as)- There 
is slightly more spread in the separate W + and W~ total cross sections. The W + /W~ 
ratio is sensitive to the u/d ratio, but insensitive to NNLO corrections (and as)- We will 
discuss W and Z production in more detail in section 6. 

4.2 gg — )■ H production 

We calculate the total cross section for Standard Model Higgs boson production (without 
decay) from gluon-gluon fusion (gg — >■ H) via a top-quark loop. The mt dependence is 
retained only at leading-order (LO), with mt = 171.3 GeV (PDG 2009 best value), and 
the higher-order corrections are calculated in the limit of an infinite top-quark mass, with 
NNLO corrections from ref. [66]. Again we use private code (originally by J. Stirling), with 
the NLO corrections taken from ref. [67] and agreement with MCFM was checked to within 
the numerical integration error (of at most 0.1%). The NNLO calculation is improved 
with respect to the numbers quoted in previous MSTW papers [48-50, 52] (which used 
mt = 175 GeV), where only terms up to (1 — z) 1 (where z = M^f-y/s) were included in the 
expansion of the "hard" pieces of the NNLO partonic cross section around the kinematic 
point z = 1 in powers of (1 — z) [66]. The improved calculation gives a 7 TeV LHC cross 
section for Mg = 120 GeV which is 1.1% larger than the expansion only up to (1 — z) 1 
and is in agreement with the exact NNLO calculation [68] to much less than 0.1%. We 
do not include the small bottom-quark loop contributions to the gg —> H cross section. 
We show the gg — > H total cross sections, plotted as a function of as(M|), in figures 11 
and 12, for Higgs boson masses (M H ) of (a,b) 120 GeV, (c,d) 180 GeV and (e,f) 240 GeV, 
with (a,c,e) 68% C.L. uncertainties on the left and (b,d,f) 90% C.L. uncertainties on the 
right. The gg — >■ H cross sections start at 0(a s ) at LO, with anomalously large higher- 
order corrections, therefore they are directly sensitive to the value of as(M%). Moreover, 
there is a known correlation between the value of as and the gluon distribution, which 
additionally affects the gg —> H cross sections. There is reasonable agreement between 
the global fits at NLO, although CTEQ6.6/CT10 give somewhat lower cross sections than 
MSTW08 /NNPDF2 . 1 , even after accounting for the slightly lower default as(Mf) value. 
This behaviour reflects the differences in the gg luminosity seen in figure 2(c,d), i.e. com- 
pared to MSTW08 and NNPDF2.1 the CTEQ6.6/CT10 gg luminosity is larger at low s 
from the positive-definite input gluon constraint and larger at high s from fitting Tevatron 
Run I jet data, leading to a smaller gg luminosity at the relevant intermediate \f§ ~ Mh 
values from the momentum sum rule. Note from figure 12 that the NNLO corrections are 
huge, and therefore essential, and are only slightly reduced by taking different as(M^) 
values at different perturbative orders. At NNLO, there is a significant discrepancy be- 
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H total cross sections, plotted as a function of as(Aff), at NLO. 



tween the MSTW08 prediction and the ABKM09 prediction, even accounting for the very 
different as(M%) values. The dependence of the gg — > H cross sections at the LHC (and 
Tevatron) on the choice of PDFs (and as) will be discussed in much more detail in ref. [29], 
where we argue that the Tevatron jet data provide an important discriminator between the 
different PDF sets and strongly favour MSTW08 over ABKM09. 
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Figure 12. — > total cross sections, plotted as a function of as(M|), at NNLO. 



4.3 it production 

We calculate ti production (without decay) for a top-quark pole mass rat = 171.3 GeV 
(PDG 2009 best value). At NLO we use private code (originally by J. Stirling) with NLO 
corrections from ref. [69], again checked to agree with MCFM to within the numerical in- 
tegration error (of around 0.1%). The NNLO calculation of the total cross section for ti 
production is still in progress (see, for example, ref. [70] for the status), although various 
approximations based on threshold resummation are available, in particular in the HATHOR 
public code [27]. We will use the HATHOR [27] code with the default settings for an approx- 
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Figure 14. tt total cross sections, plotted as a function of as(M|), at NNLO, for m t = 171.3 GeV. 



imate "NNLO" calculation, although we will make no attempt to quantify the theoretical 
uncertainty (other than from PDFs and as)', see also the discussion in ref. [71]. There 
are very slight differences at NLO because hathor [27] uses a parameterisation [72] of the 
exact analytic result [73] instead of the original numerical result [69] , giving cross sections 
around 0.03% larger (i.e. an insignificant amount) when both results are evaluated with 
high integration precision. 

We show the ti total cross sections, plotted as a function of ag(M|), in figures 13 and 
14, with (a) 68% C.L. uncertainties on the left and (b) 90% C.L. uncertainties on the right. 
More than 80% of the NLO ti cross section comes from the gg channel for the LHC with 
y/s = 7 TeV, rising to almost 90% at y/s = 14 TeV, compared to less than 15% at the 
Tevatron (y/s = 1.96 TeV). The significant difference in the initial parton composition for 
ti production is due partly to the lower Tevatron energy (pp collisions at y/s = 1.96 TeV 
would give around 50% of the ti cross section from the gg channel) , but mainly due to the 
valence- valence nature of the qq —> ti channel in pp collisions. The partonic subprocess is 
0(a s ) at LO. There is therefore a strong dependence on both the gluon distribution (at 
x ~ 2nit/y/s = 0.05) and as- The approximate NNLO corrections seem to be reasonably 
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small, especially when taking different as(M|) values at different perturbative orders, but 
there are currently many "NNLO" choices and only one possibility is shown in figure 14. A 
more complete study of the theoretical uncertainties in the approximate NNLO calculation 
is clearly important, but it is beyond the scope of this paper. Additionally, the predicted tt 
cross section has a fairly strong dependence on the assumed top-quark mass mj, such that 
comparison of the measured cross section with theory predictions even allows an extraction 
of nit [72]. As some indication of the nit dependence, moving from mt = 171.3 GeV (PDG 
2009 best value) to mt = 173.3 GeV [74] decreases the tt cross section at the 7 TeV LHC 
by about 10 pb (or 6%) at both NLO and NNLO with MSTW08 PDFs. Moreover, the 
measured tt cross section, particularly at the Tevatron, is commonly used to constrain new 
physics contributions, therefore it is questionable whether it should be used directly as a 
PDF constraint. Rather, we would hope that the gluon distribution (and as) would be 
sufficiently constrained by other data sets, such as Tevatron jet production, that the tt cross 
section is a prediction rather than a direct PDF constraint. Bearing these caveats in mind, 
the preliminary ATLAS combined measurement with 35 pb" 1 of data is 180 ± 18 pb [8] 
and the preliminary CMS combined measurement with 36 pb" 1 of data is 158 ± 19 pb [9], 
so the approximate NNLO prediction using MSTW08 shown in figure 14 is consistent 
with both ATLAS and CMS, while the central value using ABKM09 is slightly less than 
1-0- below CMS and more than 2-a below ATLAS. The discrepancy for ABKM09 would 
increase further if using a more up-to-date value of mt = 173.3 GeV [74] rather than 
mt = 171.3 GeV (PDG 2009 best value), which would reduce all theory predictions by 
around 6%. We repeat that in another paper [29] we show that the Tevatron data on jet 
production also show a definite preference for MSTW08 over ABKM09. 

5 Discussion of PDF4LHC Interim Recommendations 

The PDF4LHC Interim Recommendations [20] were originally formulated by the PDF4LHC 
Steering Committee in June 2010 [18] in response to a request from the LHC Higgs Cross 
Section Working Group [17]. The NLO and NNLO prescriptions for PDF+a5 uncertainties 
at 68% C.L. are summarised below [20]. 

NLO : "For the calculation of uncertainties at the LHC, use the envelope provided by the 
central values and PDF+a s errors from the MSTW08, CTEQ6.6 and NNPDF2.0 
PDFs." "As a central value, use the midpoint of this envelope." 

NNLO : "As a central value, use the MSTW08 prediction." Rescale the MSTW08 NNLO 
uncertainty by "the factor obtained by dividing the full uncertainty obtained from the 
envelope of MSTW, CTEQ and NNPDF results at NLO by the MSTW uncertainty 
at NLO". This rescaling factor was found to be ~ 2 for gg — > H at the LHC. 

These recommendations were motivated by the fact that the three global fits (MSTW, 
CTEQ and NNPDF) are generally in reasonable agreement, as seen in this paper, but the 
deviations between predictions are sometimes as large as their uncertainties for certain 
observables, not always for well- understood reasons. Other fits to more limited data sets, 
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particularly those not including Tevatron jet data, generally have larger deviations, again 
as clearly seen from the results of this paper. At NNLO, only the MSTW group currently 
provides NNLO PDFs, so the rescaling factor was introduced to expand the MSTW uncer- 
tainty in a similar way as at NLO, given that (experimental) PDF uncertainties are largely 
independent of perturbative order, and the observed differences between the three global 
fits (MSTW, CTEQ and NNPDF) at NLO would be expected to be similar also at NNLO. 

The PDF4LHC Interim Recommendations [20] summarised above are already some- 
what dated in that CTEQ6.6 [38] has now been replaced by CT10 [24] and NNPDF2.0 [46] 
has been replaced by NNPDF2.1 [25], in both cases correcting known defects. In particular, 
CTEQ6.6 includes only Tevatron Run I jet data, which are known to prefer a larger high-x 
gluon distribution than the more reliable Run II data [37], while NNPDF2.0 uses the in- 
adequate ZM-VFNS for DIS structure functions leading to systematically low predictions 
for W and Z cross sections at the LHC [19]. Clearly a minimal update to the prescription 
would retain the original form, but with these replacements, i.e. CTEQ6.6 — > CT10 and 
NNPDF2.0 NNPDF2.1. 

However, the key question is whether we need an official recommendation at all. The 
PDF4LHC recommendations [18, 20] were formulated by only the Steering Committee and 
should not be taken to represent a consensus decision of the PDF fitting community, as 
might be assumed. Moreover, it is difficult to see that any sensible consensus can ever be 
reached, given that many proponents of a given prescription have a vested interest that 
their own particular PDF set is used as widely as possible, thereby encouraging the use of 
"envelope" prescriptions to satisfy the majority of those involved in the decision-making 
process. The importance of PDF uncertainties is process- dependent: it is overkill to use 
a complicated prescription universally, for example, when other theoretical uncertainties 
such as scale dependence are overwhelmingly dominant. For some processes, uncertainties 
due to as are not important (for example, in the — > i^v charge asymmetry), and 
therefore the additional effort required to compute a "PDF+as'" uncertainty rather than 
a "PDF only" uncertainty is not worthwhile. For other processes, it may be necessary to 
include additional PDF-related uncertainties beyond the usual PDF+ag prescription. This 
is the case, for example, for processes initiated by bottom quarks such as i-channel single- 
top production, qb — > q't, or (supersymmetric) Higgs boson production from bottom-quark 
fusion, bb — > where the additional uncertainty due to the choice of the bottom-quark 
mass mfe can be included for the MSTW08 PDFs according to the prescription presented 
in ref. [48]; see also the NNPDF2.1 studies [25]. 

If there are large differences in predictions obtained with PDFs from different groups, 
these should be understood and results should be shown with multiple predictions, rather 
than simply taking the envelope. On the other hand, if there are only small differences 
between the predictions and their uncertainties, the envelope method will give a similar 
result to the uncertainty of any one group, so it is redundant. Moreover, the frequent nature 
of PDF releases means that any official recommendation will soon be outdated, as already 
seen with the use of CTEQ6.6 and NNPDF2.0. Indeed the updated CT10 and NNPDF2.1 
sets generally give predictions somewhat closer to MSTW08 than their predecessors. One 
possible suggestion for a simpler recommendation is to use only MSTW08 uncertainties 
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at 90% C.L. rather than the "envelope of 68%" currently proposed. This is a common 
prescription used, for example, in the 2010 Tevatron Higgs exclusion analysis [75], but not 
the most recent one [76] where the PDF4LHC recommendations have been adopted. It 
can be applied consistently at NLO and NNLO and gives similar results to the PDF4LHC 
recipe; see figure 15 where we show the gg — > H total cross section at the LHC versus 
the Higgs mass Mjj, at both NLO and NNLO, and with PDF+as uncertainties at both 
68% and 90% C.L. This simpler recommendation could be reconsidered when CTEQ and 
NNPDF sets appear at NNLO, but these are expected to follow the same trend as at NLO, 
i.e. to be in reasonably good agreement with MSTW when using 90% C.L. uncertainties. 

Despite the criticism above, the PDF4LHC prescription is not unreasonable when 
a single conservative estimate of PDF uncertainties is necessary (although it would be 
worrying if the result depended strongly on the particular recipe). This is the case, for 
example, in searches for new physics, in the calculation of acceptance corrections, and 
in the measurement of the W boson mass at hadron colliders [77], and was the type 
of situation for which the PDF4LHC recommendations were formulated. However, it is 
inappropriate for Standard Model measurements sensitive to PDFs, such as W, Z and jet 
cross sections, where the PDF4LHC recipe has commonly been applied by CMS [2, 4, 78, 
79], but not by ATLAS [1, 3, 80, 81]. Showing only the envelope masks rather than exposes 
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potential discrepancies between PDF sets, contrary to a primary motivation for these types 
of measurements as a potential PDF constraint. Indeed the original recommendations noted 
that [20]: 

"...it is useful to highlight the differences between two use cases: (1) cross 
sections which have not yet been measured (such as, for example, Higgs pro- 
duction) and (2) comparisons to existing cross sections. For the latter, the 
most useful comparisons should be to the predictions using individual PDFs 
(and their uncertainty bands). Such cross sections have the potential, for ex- 
ample, to provide information useful for modification of those PDFs." 

The PDF4LHC recipe was formulated only for case (1). Due to the intrinsic complica- 
tion such that it is difficult for non-experts to understand, let alone implement unambigu- 
ously in practice, the PDF4LHC recommmendations have already frequently been modified 
or misinterpreted, but still labelled as following the "PDF4LHC" recipe. For example, the- 
ory comparisons with CMS data on dijet production [78] and inclusive jet production [79] 
use CT10 rather than CTEQ6.6, construct the envelope and midpoint without as uncer- 
tainties, then vary as(M^) by ±0.002 using only CT10 (rather than the recommended 
±0.0012 for 1-a). Theory comparisons with the CMS data on isolated prompt photon pro- 
duction [82] obtained the uncertainty by constructing the envelope using only the central 
values of CT10, MSTW08 and NNPDF2.0, which is arguably a worse prescription than 
using only the PDF uncertainties of any one group. Envelope prescriptions for the uncer- 
tainty also have the disadvantage that PDF correlations between cross sections, such as 
those defined in section 2.2, cannot be easily calculated. 

We conclude that rather than blindly follow the PDF4LHC recommmendations, or 
some ad-hoc modification, an appropriate treatment of PDF uncertainties should be chosen 
depending on the process under consideration and the PDF sets available at the time. Of 
course, in some cases, but by no means all, constructing an envelope and midpoint may well 
be a sensible prescription, but this need not be applied universally without justification. 
In many other cases, it is better to instead consistently evaluate the complete uncertainty 
using the prescription provided by one (or preferably, more than one) PDF fitting group. 

6 Discussion of W and Z production at the LHC 

In table 3 we give percentage uncertainties at 68% and 90% C.L. of the NNLO predictions 
for W ± (= W + + W~) and Z° total cross sections, and their ratio R wz , using MSTW 2008 
PDFs [37, 48, 52]. In table 4 we show the same results for the W + and W~ total cross sec- 
tions, and their ratio R±. The asymmetric PDF [37], PDF±a 5 [52] and PDF±a s ±m Cjb [48] 
uncertainties are calculated using the prescriptions given in the respective papers. The 
uncertainties due to the choice of heavy-quark masses (m c and m b ) are estimated, for 
fixed as(M^), by varying the charm-quark mass in the range m c = 1.40 ± 0.15 GeV at 
68% C.L. and m c = 1.40 ±0.25 GeV at 90% C.L., and the bottom-quark mass in the range 
m b = 4.75 ± 0.25 GeV at 68% C.L. and m b = 4.75 ± 0.50 GeV at 90% C.L., then adding 
the variation in quadrature with the PDF±as uncertainty. The uncertainties coming from 
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+0.1 
-0.1 


+0.2 
-0.2 


+0.1 
-0.1 


+0.2 
-0.3 


+0.0 
-0.0 


+0.0 
-0.0 


PDF+Q 5 +m Cjfe 


+2.9 
-2.3 


+5.2 
-4.6 


+2.7 
-2.2 


+4.7 
-4.6 


+ 1.1 
-0.8 


+2.4 
-1.5 


Scales (fJ,R, hf) 
PDF+as+m C)b +fi RtF 


+0.5 
-0.8 
+2.9 
-2.5 


+0.9 
-1.3 
+5.3 
-4.8 


+0.6 
-0.8 
+2.7 
-2.3 


+ 1.0 
-1.3 
+4.8 
-4.8 


+0.1 
-0.1 
+1.1 
-0.8 


+0.2 
-0.1 
+2.4 
-1.6 



Table 4. Percentage uncertainties at 68% and 90% C.L. of the NNLO predictions for W + and W 
total cross sections, and their ratio R ±1 using MSTW 2008 PDFs [37, 48, 52]. 

neglected higher-order corrections beyond NNLO are estimated at 68% C.L. by allow- 
ing the renormalisation and factorisation scales (with central values fiR = hf = Mw,z) 
to vary independently so that the ratios /-Ir/Mw,z, P^fI^w,z and /j,r/hf all lie in the 
range between 1/2 and 2, then the 90% C.L. uncertainty is obtained by multiplying by 
1.64485. Finally, this uncertainty from scale variation is added in quadrature with the 
PDF+as+m C: b uncertainty, but in all cases the final uncertainty is not much more than 
the PDF+cks uncertainty used elsewhere in this paper. 

In figure 16 we show (= W + + W~) versus Z° total cross sections and in figure 17 
we show W + versus W~ total cross sections. In figures 16(a,b) and 17(a,b) the inner 
error bars represent "PDF only" uncertainties where possible, i.e. with the exception of 
ABKM09 and GJR08/JR09, and the outer error bars represent "PDF+as" uncertainties. 
The ellipses in figures 16(c,d) and 17(c,d) are drawn using eqs. (2.8) and (2.9) with the ap- 
propriate correlation cosine, where the PDF+as uncertainties have been symmetrised when 
necessary. For MSTW08, CTEQ6.6, CT10, CT10W, HERAPDF1.0 and HERAPDF1.5, 
the correlation cosine defining the error ellipse is calculated only from the eigenvector PDF 
sets for the central as value, since a definite prescription is not clear for the calculation 
of the correlation cosine in the presence of as uncertainties (and model and parameteri- 
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Figure 16. W ± (= W + + W~) versus Z° total cross sections at (a) NLO and (b) NNLO, then 
(c,d) the same plots with ellipses accounting for PDF correlations between the two cross sections. 



sation uncertainties for HERAPDF1.0 and HERAPDF1.5). We draw dotted lines where 
the ratio of cross sections, W^/Z or W + /W~ , is constant. It can be seen from the plots 
that the uncertainty in the ratio of cross sections is closely related to the correlation of the 
ellipse, also seen explicitly from eq. (2.11). In figures 16(b,d) and 17(b,d) we also show the 
NLO predictions (open symbols and dashed lines) together with the corresponding NNLO 
predictions (closed symbols and solid lines) to explicitly demonstrate how the size of the 
NNLO corrections depends on the PDF choice. 

Note that the Z cross sections presented in this paper are calculated in a zero-width 
approximation and neglecting the virtual photon (7*) contribution, whereas the LHC ex- 
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Figure 17. W + versus W total cross sections at (a) NLO and (b) NNLO, then (c,d) the same 
plots with ellipses accounting for PDF correlations between the two cross sections. 



perimental measurements include the 7* and are defined in a finite range of the invariant 
mass of the lepton pair from the Z/-f* — > decay. We can calculate correction factors 
using the vrap program [83] with MSTW08 PDFs. The results are the same at NLO and 
NNLO to three decimal places: 

ATLAS cuts [1, 3]: cr(Z/ 7 *, 66 < M u < 116 GeV) / cr(Z-only, M u = M z ) = 0.991, 
CMS cuts [2, 4]: cr(Z/ 7 *, 60 < M u < 120 GeV) / cr(Z-only, M u = M z ) = 1.006, 

i.e. there is a 1.5% difference in the predictions for the ATLAS and CMS kinematic cuts. 
Plots comparing to the preliminary ATLAS [3] and CMS [4] data can be found in ref. [23]. 
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The (dominant) luminosity uncertainty of 3-4% [3, 4] is comparable with the PDF spread 
and is much improved compared to the previous luminosity uncertainty of 11% for the 
published cross sections [1, 2], but given the preliminary nature we choose not to directly 
compare to data in this paper, particularly due to reasons associated with acceptance 
corrections which we explain below. 

The W ± and Z° total cross sections in figure 16 are highly correlated, which can 
be understood by considering the dominant partonic contributions arising from u and d 
quarks, i.e. 

cr w ++cr w - ^ u(x\) d{x 2 ) + d(x{) u(x 2 ) ^ u{x{) + d(x\) 

o^o 0.29u(ki)«(x 2 ) + 0.37d(xi)d(x 2 ) ~ 0.29 u{x x ) + 0.37 d( Xl ) ' 

where we have neglected the contributions with q o q, assuming that q{x\)q{x 2 ) dominates 
over q(xi)q(x 2 ), and the numerical values in the denominator are the appropriate sums of 
the squares of the vector and axial-vector couplings. In the last simplification in eq. (6.1) 
we have assumed that u(x 2 ) ~ d(x 2 ). Here, the momentum fractions are 

xi = (Mw/y/s) exp(+y w ) > x 2 = (M w /y/a) exp(-yw), (6.2) 
x\ = (M z /\fs) exp(+y z ) > x 2 = {M z /\fs) exp(-y z ), (6.3) 

and yw or yz should be interpreted as some "average" rapidity appropriate for the (rapidity- 
integrated) total cross section. The combination of u- and d-quark contributions is very 
similar (in numerical prefactors, in x values, and in Q 2 = M^ z values) in both the nu- 
merator and denominator of eq. (6.1), therefore the PDF dependence almost cancels in the 
W^/Z ratio. The W + and W~ total cross sections in figure 17 are much less correlated, 
since 

&W+ ^ u(xi)d(x 2 ) ^ u(xi) 
o w - d(xi)u(x 2 ) d(xi)' 

and therefore the W + /W~ cross-section ratio is a sensitive probe of the u/d ratio and 
there is more variation between different PDF sets. The preliminary CMS measurement of 
the W + /W~ total cross-section ratio with 36 pb _1 of integrated luminosity (the full 2010 
data set), combining the electron and muon channels, is [4] 

= 1.421 ± 0.006(stat.) ± 0.014(syst.) ± 0.030(th.). (6.5) 

<7W- 

Comparing to the predictions shown in figure 17, or more easily those in figure 10, we see 
that the preliminary CMS measurement is in good agreement with MSTW08 and is more 
than 2-cj below the central value of ABKM09. 

The (dominant) theory uncertainty in the preliminary CMS W + /W~ total cross- 
section ratio, eq. (6.5), comes from extrapolating the measured (fiducial) cross section 
over the whole phase space, where the central value of the acceptance is obtained using the 
powheg Monte Carlo event generator with CT10 PDFs (with MSTW08 and NNPDF2.0 
PDFs used to evaluate the uncertainty). The corresponding ATLAS measurement of the 
W + /W~ ratio with 35 pb _1 of integrated luminosity has not yet been explicitly presented. 
The acceptance corrections for the central values of the individual preliminary cross-section 
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Figure 18. ATLAS W -> tv charge asymmetry [85] compared to MSTW08 NLO PDFs with 90% 
C.L. uncertainty band, (a) Comparison also to MRST LO* (used as the central value for the ATLAS 
acceptance calculation [3]), CTEQ6.6 and HERAPDF1.0 (used to evaluate the uncertainty [3]). 
(b) Comparison to various other LO, LO* and NLO PDFs from MRST/MSTW and CTEQ/CT. 



measurements [3] were obtained using pythia with MRST LO* PDFs [84], and seem to 
give an anomalously low value of acceptance ratios leading to a high value of the central 
W + /W~ total cross-section ratio of 1.51 [3], whereas alternative acceptance calculations 
would move the central ATLAS W + /W~ total cross-section ratio much closer to the CMS 
result [23]. Moreover, the MRST LO* PDFs [84] give a poor description of the ATLAS 
data [85] on the W — > Iv charge asymmetry inside the acceptance, worse than any other 
PDF set (see figure 18). The W —> iv charge asymmetry is defined differentially as a 
function of the pseudorapidity rjg of the charged-lepton, i.e. 

Mm) = d^+)/d^ + d^-)AV (6 - 6) 

Here, common NLO K- factors have been applied for all curves in figure 18, but these have 
only a very small effect [23] . The LO* notation denotes a modified LO fit with use of a NLO 
strong coupling and relaxation of the momentum sum rule [84] , the latter feature perhaps 
leading to less constraint on the relevant u/d ratio at x ~ Mw/^/s ~ 0.01. The analogous 
CT09MC2 [86] prediction also differs substantially from the corresponding CTEQ6L1 and 
CTEQ6.6 predictions; see figure 18(b). It is clear that the current data precision requires 
better acceptance calculations for the W and Z total cross sections and/or data to theory 
comparisons at the level of the fiducial cross section [23]. 

Fiducial cross sections have been quoted by ATLAS [1, 3] and CMS [2, 4], but not for 
the cross-section ratios. The LHCb experiment only quote fiducial W and Z cross sections 
(and ratios) [5] given that extrapolation to total cross sections would not be meaningful due 
to the limited LHCb acceptance. "Fully exclusive" calculations at NNLO, allowing direct 
computation of the fiducial cross section within leptonic cuts, are available in the form of 
the public fewz [87] and dynnlo [88] codes. However, data-to-theory comparisons using 
these tools at the level of detail performed for total cross sections in this paper would be 
arduous (see refs. [89-91] for some early NNLO fewz studies with leptonic cuts), due to the 
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difficulty 5 in achieving sufficient numerical integration precision in a reasonable amount of 
CPU time, and the need for multiple independent runs to evaluate PDF uncertainties. The 
situation has improved recently with the release of the updated fewz 2.0 [92], with refined 
integration routines and automatic calculation of PDF uncertainties in a single run, but so 
far only for Z/j* production. After the similarly-improved fewz 2.0 for production is 
released, then detailed NNLO comparisons at the level of the fiducial cross section will be 
more practical, although the leptonic cuts are different for ATLAS and CMS (and LHCb), 
and for the electron and muon channels, so the data-to-theory comparisons would need to 
be performed separately for each experimental measurement. 

The large spread of predictions for the W + /W~ ratio, and the related W — > Iv charge 
asymmetry, between different PDF sets is of some concern. We recalculated the W + /W~ 
total cross-section ratios with only u and d quark (and antiquark) distributions to confirm 
that the same trend is seen. To examine sources of tension we performed a variant of the 
MSTW08 NLO fit with a pseudo-data point added corresponding to the CTEQ6.6 value 
of the W + /W~ total cross-section ratio at the 7 TeV LHC (see figure 9) with a small 
uncertainty. The data sets for which the fit quality most deteriorated were the Tevatron 
Run II W — > Iv asymmetry [93, 94], the NMC fin/up ratio [95] and the E866/NuSea pd/pp 
Drell-Yan ratio [96]. The CTEQ6.6/CT10 input parameterisation, and that of some other 
groups such as HERAPDF, assumes the same small- x behaviour of the u and d valence- 
quark distributions, whereas the MSTW08 NLO fit has independent small-x powers, finding 
that at the input scale Qq = 1 GeV 2 [37]: 

xu v (x,Ql) oc x°- 29 -o X oi, xd v (x,Ql) oc x°' 97 ±o.ii. (6.7) 

The additional freedom in parameterisation may provide some explanation of the difference 
between MSTW08 and CTEQ6.6 at central pseudorapidities seen in figure 18 and for 
the W + /W~ total cross-section ratio. On the other hand, the preference of the ATLAS 
W — > Iv charge asymmetry data [85], and also of the CMS W —> iv charge asymmetry 
data [97], seems to be for CTEQ6.6/CT10 rather than MSTW08. This implies some tension 
between the LHC W —> Iv charge asymmetry data [85, 97] and the data already included 
in the MSTW08 fit. Other tensions have been observed with the precise Tevatron data on 
the W — > Iv charge asymmetry, and partially resolved by more flexible nuclear corrections 
for deuteron structure functions [50]. Further attempts to resolve these tensions will be 
necessary for any future update of the MSTW08 fit and this will not be a straightforward 
endeavour. Other than the LHC W — )• Iv charge asymmetry data [5, 85, 97], further PDF 
constraints from electroweak data will come from the Z/^f* rapidity distribution, already 
presented in preliminary form by CMS [98] and LHCb [5]. 

7 Conclusions 

We presented results on total cross sections versus as(M^) for the LHC (y/s = 7 TeV) 
benchmark processes: W ± , Z°, ti, gg ->■ H (M H = 120,180,240 GeV). This study was 

This technical difficulty meant that it was not possible to compute accurate NNLO corrections using 
fewz for the Tevatron W -> iv charge asymmetry data included in the MSTW 2008 NNLO fit [37]. 
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an extension and update of the results presented in ref. [13]. All plots (and more) are 
available from a webpage [28], which may be updated in future as more results become 
available. Comparing predictions at the same as(M^) value allows a fairer comparison, 
but there are still sizeable differences between results from some PDF groups, particularly 
those groups using more limited data sets and/or restrictive input PDF parameterisations. 
Looking at partonic luminosities can often give a more global picture of PDF differences 
compared to total cross sections. Similar differences are apparent by looking at rapidity 
distributions (see figures 11-13 of ref. [13]). We did not, in general, attempt to evaluate 
theory uncertainties other than those arising from PDFs and as, which roughly decouple 
from other considerations such as scale dependence, therefore the results presented here 
will be valid also for more complete studies. 

There is now reasonably good agreement between the NLO global fits from MSTW08 [37], 
CT10 [24] and NNPDF2.1 [25], all using variants of a GM-VFNS to treat DIS structure 
functions. There is more variation with other PDF sets, particularly at NNLO where 
MSTW08 is the only public global PDF fit available. The tt and gg —> H cross sections are 
sensitive to the gg luminosity and as(M§), but we can use Tevatron jet data to discrim- 
inate [29]. Fits not including Tevatron jet data should therefore be used with caution (if 
at all). We discussed the recent PDF4LHC Interim Recommendations [20]. A proposal for 
a simpler recipe is to use only MSTW08 with PDF+ag uncertainties at 90% C.L., which 
gives a common prescription at NLO and NNLO, with similar results to the "envelope at 
68% C.L." prescription suggested by the PDF4LHC Steering Committee [20]. The 2010 
Tevatron Higgs exclusion analysis [75], using gg — > H cross sections calculated with only 
MSTW08 at 90% C.L., is robust using exactly this prescription, although the more recent 
analysis [76] instead uses the PDF4LHC recipe which gives similar results; see ref. [29] for 
more details. 

We discussed the comparison of data and theory for W and Z production at the 
LHC. The ATLAS and CMS data [3, 4] are now reaching a level of precision where the 
theoretical uncertainty on the total cross sections due to the acceptance corrections is 
sizeable, particularly for cross-section ratios, and comparison at the level of the fiducial 
cross section would be desirable. Nevertheless, the total cross section remains a convenient 
"pseudo-observable" for examining differences due to PDFs, even if it does not allow a 
direct comparison to be made between data and theory. The more differential LHC W —> Iv 
charge asymmetry data [5, 85, 97], defined only within fiducial cuts, will be addressed in 
more detail in a future study, together with data on the Z/j* rapidity distribution [5, 98]. 
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